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ABSTRACT 

A theory is presented for computing the far field beam 
patterns from distributed random noise sources. The theore- 
tical model utilizes the Green's Function for the wave equa- 
tion and the space-time autocorrelation function for deter- 
mining the radiation from a randomly vibrating area. The 
actual far field beam pattern of a horn speaker in an anechoic 
chamber was obtained, and also near field measurements were 
taken to obtain the correlation distance and the mean square 
of the particle velocity using the autocorrelation function. 
Finally a computer program was written to evaluate the inte- 
gral wave equation by numerical methods. It was found that 
the critical parameters in the mathematical model were the 
correlation distance and the frequency limits of integration. 
Small variations in the correlation distance modified greatly 
the width of the predicted beam pattern, while changes in the 
limits of integration had a moderate effect. The Frequency 
Spectrum was obtained in the anechoic chamber and it was used 
to determine the limits of integration of the integral solu- 
tion for the intensity field. 

The predicted beam patterns were in satisfactory agreement 
with the measured far field beam patterns. The differences 
between them are probably due to inaccurates in determining 
the limits of integration and the correlation distance. 
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I. INTRODUCTION 



It has long been recognized that the best means of locat- 
ing a deeply submerged submarine is by detecting the sound 
that it radiates into the water; also all surface ships are 
subject to attack by passive homing torpedoes and acoustic 
mines, both of which operate using the radiated noise of the 
ship. These facts make the far field radiated signatures a 
very important factor to determine detectability characteris- 
tics and capacity of survival in a hostile environment. 

Free field measurement in the far field involves serious 
difficulties like ambient noise, multipath, convergence zone, 
attenuation, among others, which makes impracticable such 
kind of measurement. Instead a prediction of the Far Field 
using Near Field data is suggested. Certain general princi- 
ples are repeatedly invoked in prediction. The first is that 
the vibratory powers of several sources is additive. This is 
based upon the assumption that the signals from any two 
sources are statistically uncorrelated, an assumption which 
is almost invariably found correct. Signal contributions from 
non-overlapping frequency bands are inherently uncorrelated; 
hence the total power over all frequencies is found by summing 
contributions from various frequency bands . 

Several works have been done about prediction of Far 
Field beam pattern from Near Field measurements for discrete 
frequencies (Refs. 1,2, 3, 4). In general these methods deal 
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with the amplitude and phase information taken at the Near 
Field and the final result shows that the Fourier Transform 
of those measurements lead to the beam pattern. Consider 
the random noise as the random summation of sinusoidal waves, 
thus a complete analysis of the frequency bands can be 
achieved using a very narrow filter which center frequency 
can be shifted to cover all the spectrum, or a bank of narrow 
filters can be used. Adding the contribution of each fre- 
quency it should be possible to obtain the final beam pattern, 
but either procedure involves a tedious work, specially for 
broad band signals. There also remains the question of cor- 
relation between different bands. 

The more general problem of predicting the Far Field noise 
beam pattern from a random noise source is of great practical 
importance. The theory of the noise field for systems that 
satisfy the wave equation has been developed by some inves- 
tigators working in the fields of acoustics (Ref. 5), and 
optics (Ref. 6). Different methods had been suggested in or- 
der to generate a prediction in the far field, one of them 
consists of using correlation functions and the Green's func- 
tion for the wave equation, for general time variations. This 
technique is used by Morse and Ingard (Ref. 7) in their book 
"Theoretical Acoustic," for describing non-periodic vibra- 
tions or random processes. 

This paper discusses one possible method of Near Field 
analysis using the spatial correlation function and the 
Green’s function which enables differential equations and 
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boundary conditions to be combined into an integral equation 
for the wave motion; the solution of the integral equation 
provides the intensity in the Far Field. This method could 
practically be used in collecting submarine radiated acoustic 
signatures . 
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II. GENERAL THEORY 



The analysis of the disturbed region has to be done us- 
ing a different approach than for discrete frequencies and 
one of the most common becomes the correlation technique. 

It is well known what the relationship is between the auto- 
correlation function of random signal and the total average 
power and the D.C. power contained in the signal. In a 
sense the autocorrelation function is a measure of time 
variation, spatial variation and statistical dependence in 
the region of inhomogeneities. 

The study of the statistics of noise gives the possi- 
bility of finding the solution of the problem using the cor 
relation function instead of a frequency decomposition and 
individual analysis of the frequency components in the ran- 
dom signal. In the process any phase information is lost, 
but the important point becomes the relation between the 
autocorrelation function and the field intensity, which is 
phase independent, and the angular dependence with respect 
to the acoustic axis of the source. 

It is helpful to begin with a non periodical vibration 
problem, with several degrees of freedom, each degree of 
freedom has its own equation of motion and the general equa 
tion becomes: 

d^X ^ 

— ^ = F^(r) + F* (r,t) + D„(v,t) (2.1) 

n n n n ’ 

at 

n = 1, 2 , . . .N 
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where F (r) is the internal force, F' (r,t) is the exter- 
n ’ n 

nal force, D (v,t) is the frictional force, and M is the 
n ’ n 

effective mass of the system. 

Plotting the behavior of one component of x and its 
associated velocity on the Nth phase plane, it can be shown 
that the contour is not a closed- loop. Usually it will be a 
complicated non periodic curve confined within some limits 
given by the energy of the system but never repeating itself. 
Of course the N gegree of freedom implies N different frequen- 
cies in such a way that the combination is not periodic. The 
displacement component will not be a periodic function of 
time, so the Fourier Analysis can not be applied. The average 
value of the displacement is also time dependent for short 
periods of time, but if the time is long enough the fluctua- 
tion can be evened out, and the resultant average becomes a 
constant; again it is necessary to use the statistical pro- 
perties of the autocorrelation function, which is defined as 
follows : 

R(t) = Lim [i x(t) • x(t+T) dt] (2.2) 

-hi 

T-><» 

The value of R(0) is equal to the average power of the 
signal over the interval, and this will be the maximum value 
of R(t) . The size of "t" is very important since the 
rapidity with which R(x) drops to zero as t is increased 
is a measure of the degree of randomness of the oscillation 
of the function. In the case of the spatial autocorrelation 
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function, the correlation distance, (d), follows the same 
behavior. 

The second step will be to study the radiation character- 
istic from a plane surface with randomly vibrating area, the 
goal in this section is to try to obtain a Green’s function 
such that one or more parts of the resulting integral equa- 
tion can be eliminated, and the remaining part easy enough to 
be integrated. The solution of this problem will be approached 
in four steps : 

(1) Calculation of the far field due to a point source in 

the presence of a rigid boundary serves to outline the geometry 
and describe the role of the boundary and facilitate the choice 
of the right Green’s function for subsequent steps. 

(2) Calculation of the far field due to harmonic motion of 
parts of the surfaces itself which is assumed to have surfaces 
characteristic of a rigid boundary. 

(3) Calculation of the far field due to harmonic motion of 
parts of the surfaces which have a particular acoustic 
impedance. 

(4) Application of results from (3) to randomly vibrating area. 
1. Consider sources in the presence of an infinite bar- 
rier like the arrangement shown in Figure 2.1. If the plane 

is perfectly rigid, the boundary condition is that the par- 
ticle velocity and the gradient of the pressure at the surfaces 
is zero, thus, the effect of the boundary plane can be replaced 
by a set of images sources, symmetrically located with respect 
to the plane, both radiating into the region containing Po . 
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defined by the 



The acoustic pressure at Q 
vector r , for a simple-harmonic source is then: 



p = -ipok 




(2.3) 



where k is the wave number and represent the strength 
of the source expressed in the frequency domain. Some 
approximations have to be made in Equation (2.3), assuming 
that Q is at a distance such that the two waves combine 
to form a source that looks like a single nonsimple source 
located at the origin. The magnitude' of vectors r^ and 
r 2 is defined by: 



Important parameter since it determines the constructive or 
destructive interference. For example, if kz^ and kr^ are 



r 



1 




r - r„ cos 9 



o 



and 




r + 2z^ cos a - r^ cos (p 
o o 



Using these approximations Equation (2.3) becomes: 




lk(r - r cos 9 + z^ cos a) 
' o o 



The distance from the source to the plane is an 
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much less than 1 , the source is much smaller than a 
wavelength and the behavior of the pressure in the far field 
is very similar to that produced by a simple source of 
double strength of the original source. On the other hand, 
if is equal to a quarter wavelength the separation 

between the source and its image is one half the wavelength 
and its effect is a destructive interference. The same 
effect can be produced if instead of a point source a 
multipole source exists for which each point has its own 
image. Of course the boundary produces distortion in the 
near field which is not considered in this paper. 

2. Now, considering that the sound waves in the bounded 
medium may be produced by the motion of the plane surface 
Itself. With the plane rigid at z equal to zero, some 
portion of the plane may be moving in the z direction. 

For a simple harmonic motion u^ = u^ exp(-icot) , the 
boundary condition is that the gradient of the pressure is 
ikpcu^(x,y)exp(-iwt) , whereas, for the rest of the plane 
the partial derivative of the pressure with respect to z 
is zero. 

The time Independent part of the solution of the 
wave equation in terms of the Green's function is given by: 









G (— )] 
w r^ 
o 



dS, 



o 

(2.4) 
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since there are no sources in the medium the volume integral 
is equal to zero and the surface integral Involves p and 
the outward pointing normal gradient at z = 0 , but the 
description of the physical situation has specified only 
the normal gradient. The Green’s function should be selected 
in such a way to eliminate the second part of the surface 
integral, keeping only the specified normal gradient of the 
pressure. The Green's function satisfying this requirement 
is given by: 



G (-^) 
“ ^o 





ikr2 

e 



in this case ~ ^2 * because the source is in the surface 

boundary, thus, the bounded medium Green's function 
is twice as much as the unbounded medium Green's function 
g^ , the reason is that each moving element of area radiates 
only into the half space z > 0 , the source and its image 
have coalesced to form a source of double strength. 

Applying this concept. Equation (2.^1) becomes: 



p^(x,y,z) 



-ikpc //u^(x 
-21kpc //u^ 






( X V ) S f ^ — ) 



<^^0 ^^o 



(2.5) 



3 . A perfectly rigid boundary is unrealizable. The 
common situation is that the boundary has a particular 
acoustical impedance defined by its density and the velocity 
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of sound in the medium. Then, the boundary reacts to an 
incident pressure and every portion of the plane is moving, 
and a different method has to be applied in order to dis- 
tinguish the active regions in the plane, where energy is 
fed into the medium from the passive areas. 

For the passive areas the relationship between p 
and its normal derivative is given by: 



-u. 



ikpc 



IE 

dZ 



-^P 
pc ^ 



at z = 0 



( 2 . 6 ) 



where B is the impedance ratio between the impedance of 
the medium where sound propagates and the impedance of the 
surface boundary. For the Active Areas the normal fluid 
velocity, plus the admittance of the boundary times the 
pressure at z = 0 must differ from zero, this difference 
will be defined as the driving velocity (^^J) 



(u. 



pc 



P) 



z=0 



(-i- + 

ikpc Zz 



_L 

pc 



P) 

z=0 



u e 
0) 



-Iwt 



( 2 . 7 ) 



To simplify the surface Integral, the Green's function 
could be arranged such that the surface integral over the 
passive areas is zero. This is accomplished by a Green's 

I 

function that satisfies: 



3G 

^ ( 2 . 8 ) 
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The integral over the active area is: 



P 




(2.9) 



Using Eq. (2.8) in Eq. (2.9) 



P 



01 




( 2 . 10 ) 



The expression in parentheses is the driving velocity in 
the active region given by (2.7), thus Eq. (2.10) becomes: 



which has the same form as Equation (2.5), but this time an 
exact representation of the Green’s function as was stated 
previously is not possible. However, a good approximation 
valid when the measurement point Q is not closer than a 
half wavelength to the boundary surface and which satisfy 
Equation (2.8) is given by: 



where the variables correspond to Fig. 2.1. Although this 
form of the Green’s function is derived for the problem of 
a volume source in the presence of a non-rigid boundary, it 
may be used since it accomplishes the simplification of the 
problem in hand. [In this approximation there is an image 



-ikpo 0 




( 2 . 11 ) 




( 2 . 12 ) 
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at Pq , but the wave it radiates has an amplitude which 
depends on the angle a' which the line r2 makes with 
the z-axis normal to the boundary. The ray of acoustic 
radiation which reaches the point Q from the source point 
at Pq by reflection, strikes the reflecting surface at Q 
at an angle a’ of Incidence and reflection. The wave 
which goes directly from P^ to Q without reflection, 
given by the first term, is independent of the inclination 
of r^ ; the reflected wave, given by the second term, 
depends on the angle of reflection. The factor 3 introduced 
the variation in the reflected sound due to the boundary. In 
this way the term in parentheses represents the reflection 
coefficient . ] 

In the far field Eq. (2.12) can be simplified using 



= r. 



G 

(0 



ikr, ikr, 

e cos g* - B -| _ e r cos a* -■ 

Mur^ ^ cos a* + 3"^ 27rr^ ^cos a* + 3"^ 



(2.13) 



From Eqs. (2.11) and (2.13) the pressure in the far 
field is given by: 



0) 






ikr. 



cos g* 
cos a' + 3 



dx^ dy^ 
o o 



(2.14) 



Observe that the integral is over the active region 
of the plane which is assumed to be confined to the circle 
of radius ”a” about the origin. 



19 



In the far field the last equation can be simplified 
neglecting terms of order a/r and smaller and using the 
approximation cos a* = cos a = z/r , 



\ cos <}) + Yq sin <j>) sin a 



1 ^ 1 



and tan (}) s — 



Thus, Equation (2.7) becomes: 



Pjr) :: 



-ikpc 



cos g 



g + cos a 



ikr -ik*W sin a 

(2.15) 



where the vector lies in the xy plane, with x^ , 

as component; and vector k has the r direction whose 
components are: 



k = k cos (j) sin a 

k = k sin (j) sin a 

k = k cos a 

z 

A close study of the integral part of Equation (2.15) 
shows that the angular distribution of the Far Field is 
proportional to the two dimensional Fourier Transform of the 
distribution of the driving velocity in the active region, 
which is defined by: 
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1 




00 




dk dk 
X y 




U^(k cos <j> sin a, k sin <|> sin a) (2.l6) 



Thus Equation (2.15) takes the form: 




^ikr 

e 



(2.17) 



The angular dependence is given by the factor , which 

measures the relative amplitude of the pressure with respect 
to the center of the active region. 



If the coordinate axes are rotated to X^.Y , such that X 

o ' o ’ o 

axis is the intersection of the rz plane, the appropriate 
transformation is given by: 

^o * ^o ^ '^o 

Under this transformation the time independent particle 
velocity becomes; 



Several properties of the Par Field can be studied. 
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-ikX sin a 
o 



U^(k cos (p sin a,k sin (j) -sin a) 




®o J' ®o 



This equation shows the one dimensional Fourier 
transformation in the direction, of the average; 









dY. 



of the particle velocity u (X^,Y ) in the Y direction. 

woo o 

It can also be observed that u^(a),X^) is a function 

of X , the distance in the direction of the measured 
o 

point Q ; it is also a function of (}> , the direction of 
the rz plane, and of the frequency w/2tt . 

If the particle velocity u (x,y,t) is an arbitrary 

Zt 

function of space and time and u (x,y) is its time dimension 

w 

Fourier transform, then the pressure field at a' point Q (r,a,(})) 
is the inverse transform of given by Equation (2.17) 

00 — ICl) ( t — — ) 

p(r,t) = / f^(aj(j)) e (2.18) 

“00 



Now, applying the rotation of axes as defined before. 
Equation (2.18) can be expressed as: 



nfr t') = -5— . cos g r 

2irr 6 + cos a ^ 



cos a 



<» -lio(t - 'I’) 



-lioe 



~ 1(^)X sin a 
'c o 

d(D / e 



dX 



^X> 

00 






(2.19) 



p(r,t) = • a' Itt [/ dX /dY u(X ,Y ,t + ^ X sin a - J)] 

* 2rrr 3 + cos a9t o o o*o* co 



P? 



This equation shows clearly some of the properties of the 
Far Field. An appropriate summation of the normal acceleration 
of the active area gives the total field at the point of 
Interest. Also the contribution to the pressure at r , at 
time t , from the line , -«> < < oo ^ depends on 

the normal acceleration on the line at time 

t+(l/c)XQ sln(a)-r/c , at time r/c - (X^/c) sln(a) earlier, 
which Is the time required for the sound to travel from the 
line to point Q . 

4. All the formulae developed so far can be used to 
compute the mean Intensity at r,a,4», produced by an area A 
of the boundary which Is vibrating with a random distribution 
of velocity. The pahtlcle velocity In the active area A Is 
a function of position as well 'as of time. In this regard. 

It becomes a random variable and a way to describe Its 
behavior Is using Its statistical properties. In this 
particular case the autocorrelation function Is used as a 
means to describe the behavior of the random source. 

Considering the particle velocity as a function of position 
and time in the active area, the autocorrelation function 
is defined by: 

H(C,n,T) = ^ /// u*(x,y,t) u^(x+C,y+n,t+T) dx dy dt 

where the integration is over the active area A , and T 
is over a long interval of time. Also the Fourier transform 
of the velocity is given by: 
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lojt-ik^x-lkyy 



dx dy dt . 






Prom this equation the space time spectrum level of 
the particle velocity is: 




AT 

6Mit^ 



±wT-ik £-ik n 

/// R(C,n,T)e ^ y dCdndT 



( 2 . 20 ) 



In the far field it is known that the intensity is 
given by the real part of the pressure times the particle 
velocity, but also the velocity is related to pressure for 
the relation u = p/pc , thus using Eq. (2.16) and Eq. (2.17) 

2 2 2 2 
|-^1^ s: ^ — — — I |U (k sin a cos k sin a sin <p) \ 

r'^ I g + cos a I ^ 



( 2 . 21 ) 



This equation is related to the autocorrelation 
function by Equation (2.20). The problem is now reduced to 
getting the form of R(C>n,T) in space and time. A 
reasonable assumption could be a Gaussian waveform where the 
value of R falls off smoothly as its arguments are Increased 
Assume that the autocorrelation is given by: 



R( 5>n>T) 




, -2 ^ 2x 1 / s2 

( C + n ) - 2 



( 2 ^. 22 ) 



2 

where is the mean square particle velocity of the active 

area averaged over time interval T ; d stands for the 
correlation distance along the plane, is the upper 
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angular frequency limit of the spectrum, which also can 
give the correlation time, and the lower limit of the 
wavelength . 

Using this value of the autocorrelation function 
and equations (2.20) and (2.21), the intensity field was 
evaluated: 



I^(w) t pc ^ 
r 



2 2 ^ 
c (0_l3 + cos al 
n* ‘ 



But this is only for a specific frequency; in order 
to get the total field intensity it would be necessary to 
consider the contribution of all the frequency components 
of the random signal. A logical way to accomplish this is 
to integrate Equation (2.23) with respect to frequency. 

Note that the integrand has a maximum when the exponent 
is zero at: 



w = 2 ^'^^ 0)^(1 + sln^ 



a) 



1/2 



GO 




/ Ij,(w)dw s 




|Uq|^ cos^ a (kj^d)2 

[3 + cos a|^ [1 + (kj^d)^sin^ 



(2.24) 



A simple analysis of this formula shows that the 
Intensity is proportional to the Impedance of the medium. 
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times the mean square particle velocity in the active area, 
times the ratio between the active area, A , and the area 
of the half sphere of radius r , times an angle distribution 
function , which depends on the maximum frequency of the 
spectrum (inverse of the correlation time), the speed of 
sound and the correlation distance, d . 

Trying to apply this theory to practical cases, two 
main problems arose. The first is that in real situations 
the frequency spectrum is band limited, usually due to the 
physical characteristics of the source. Under this 
consideration the limit of integration should have an upper 
bound (UL) and a lower limit (LL) , in this form. Equation 
(2.23) becomes: 



I (r) = 2pc ^ (27t)"^/^ 
^ r 




^tod,.2 . 2 . 1 , (0 x2 ,,, 

sin do 

oT e ^ 



(2.25) 



There is no analytical solution for this equation, so 
a numerical method of integration had to be used in order to 
get the intensity in the Far Field. 

The second problem is the difficulty in measuring the 
autocorrelation function of the randomly vibrating area, 
which is related to the particle velocity in the active area. 
Instead it can be assumed that if the area close to the 
active region is vibrating with similar characteristics, 
then the impedance ratio becomes unity and the autocorrelation 



26 



in the Near Field gives the mean square particle velocity; 
Equation (2.25) becomes: 






2po C2u) 



.3/2 d^lUpl^ cos^ g 



2 2 
c 03 ll + COS a I 
n' ' 



UL 
/ to^e 
LL 



^ojdv2 .2 1/ 0) n 2 j 

- ^— ) sin a - dw 

2 c 2 0 )^ 



( 2 . 26 ) 



This equation is angle dependent, thus the intensity 
in the Far Field depends on the angular position with respect 
to the acoustical axis. In this way it is possible to 
evaluate the beam pattern. Observe also that should 

lay inside of the limit of integration or at least equal to 
the upper limit, otherwise the process would show a 
discontinuity which is not possible because the frequencies 
contributing to the field are given by the limit of 
Integration. 
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III. EXPERIMENTAL RESULTS USING THE PLANE MODEL 



To validate the theoretical model described in Section 
II, Far Field and Near Field measurements were made on a 
horn speaker in an anechoic chamber. The actual beam pattern 
in the far field of the horn-speaker was taken and also the 
autocorrelation function in different points of the near 
field were obtained in order to get the mean square particle 
velocity in the Near-Field and be able to predict the Far 
Field Beam Pattern. 

A. FAR- FIELD MEASUREMENT 

The normalized Far Field beam pattern was obtained using 
the experimental set up shown in Figure 3.1. Since the 
source was driven by random noise the distance of the near 
field is variable. The criterion for near field is usually 
given for one particular frequency, as the frequency is made 
lower and lower the near field also decreases and eventually 
disappears for practical purposes. That is the reason for 
the presence of the high pass filter which limits the low 
frequency component to 1300 Hz in order to get an appreciable 
near field region. This is necessary for obtaining far field 
and near field measurements under the same frequency spectrum. 
The signal is then amplified and applied to the horn speaker. 
The microphone was placed in the far field at a distance of 
5o8 meters from the source, the received signal was amplified 
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Figure 3.1 Far field beam pattern measurement setup. 




and the output connected to a true RMS Voltmeter which also 
generates a DC voltage output proportional to the deflection 
of the needle for the selected scale, this voltage was ap- 
plied to the Y axis of the plotter. The X axis of the 
plotter received information from the position indicator as 
a DC voltage proportional to the angular displacement of the 
acoustic axis of the horn speaker. The normalized pressure 
amplitude versus angular position is shown in Figure 3.2. 

This was the actual Far Field Beam Pattern which had to 
be compared with the beam pattern evaluated using the mathe- 
matical model described in the previous section. 

B. NEAR-FIELD MEASUREMENT 

The intensity in the Far Field, as defined by Equation 
(2.26) depends, among other factors, upon the mean square 
particle velocity. The mean square particle velocity is a 
source quantity which was shown to be related to the auto- 
correlation function. The evaluation of the factor is 
possible by using an array of microphones over the entire ac- 
tive region and the summation of the received signal fed 
into the correlator. The microphone array was not available, 
instead a single microphone was used to take a measurement 
each 5 cm in the vertical and horizontal direction trying to 
cover the active area as shown in Figure 3.3. The autocorre- 
lation values obtained at each single point is then averaged 
and the resultant value was used to describe the mean square 
particle velocity. 
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Figure 3.2 Measured far field beam pattern. Horn Speaker. 
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Figure 3.3 Near field schematic representation of 

the position of the microphone. (View facing horn) 



The correlator was a SAICOR model SAI-43A, which was 
used in the autocorrelation mode, the function is determined 
simultaneously at 400 incremental lag points, the completed 
correlation function is displayed at one time, the function 
is provided with incremental lag or time delay value rang- 
ing from 0.2 microseconds to 1 sec resulting in total time 
delays from 80 microseconds to 400 sec. Precomputation de- 
lay of 800 lag values selected in 200 lag increments allows 
the autocorrelation function to be viewed symmetrically 
about zero or up to 800 lag values removed from zero. The 

averaging is accomplished digitally in a linear fashion with 

9 17 

fixed summations ranging from 2 to 2 in binary steps. 

The quantization is 8 bits per channel. 

The SAI-43 operates as follows: one input signal is 

applied to either channel, and internally routed to the 
second channel, the signal is sampled at the rate indicated 
by the operator-selected sample increment and then is quan- 
tized in both channels. Channel I is the channel which con- 
tains the 400 word memory. Channel II contains the one word 
memory - the present sample. At the time called "Present 
Time", the correlator performs the multiplication of the 
present value word of Channel II by each of the 400 words 
in Channel I. The product "0" times "0" becomes an entry in 
R(0), the first point of the autocorrelation function the 
product "0" times "1" becomes an entry for R(l) and finally 
"0" times "399" becomes an entry for R(399). At the end of 
this arithmetic cycle the correlator has performed one calcu- 
lation for each so-called bin of the autocorrelation function. 
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The product in R(0) is the first entry for the zero-lag auto- 
correlation value, while the product in R(399) is the first 
entry for the correlation value at a lag in Channel I equal 
to 399t. These first entries, resulting from one sample in 
Channel II are expanded to the full correlation function by 
taking repeatedly updated value of the sample in Channel II 
and correspondingly updated values in Channel I and then 
averaged. In effect the 400 word memory slides ahead one 
word and the arithmetic cycle begins again. For all t incre- 
ments that are equal or longer than 500 usee the SAI-43A is 
said to operate in real time. But for increments less than 
500 usee the figure is different. Assume an increment less 
than 50 msec, it could have completed eight cycles in 400 
msec instead of one, therefore it is eight times slower than 
the real time instrument. 

The output of the correlator can be displayed in an 
oscilloscope or in an x-y plotter, the voltage output in the 
autocorrelation mode can be determined by the following 
formula : 

r(„ 2) . 2(^ao)dBA/10 (3.^, 

where x is the voltage as seen as the output of the SAI-43A 

dBA is the dB setting of the input channel. R is the 

. 2 

correlatxon value in (volt) to be determined. 

For illustration purposes Figure 3.4 and Figure 3.5 show 
the autocorrelation function, R( 0,0,0), for the point of 
maximum intensity, (center of the speaker), and the point of 
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Figure 3*^ Autocorrelation function measured at the boundary 
of the active region 
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Figure 3.5 Autocorrelation function measured at the center of the speaker 



minimum intensity, which correspond to the corner of the 
boundary of the active region, the relative power between 
them is more than 20 dB. 

The physical arrangement in the anechoic chamber and the 
instrumentation for taking the autocorrelation function at 
each point in the near field is shown in Figure 3.6. The 
autocorrelation of the band limited noise becomes the sin x/ 

X, function, the amplitude of the function is variable depend- 
ing on the position of the microphone. As described before 
it is also possible to get the power at each point of the 
autocorrelation function. 

The criteria used for determining the active area was on 
the basis of the differences in power between the center of 
the horn up to the point which had 20 dB less power, which 
corresponds to the boundary of the active area shown in 
Figure 3.3. 

The other important measurement was the frequency spectrum 
of the received signal which is given in Figure 3.7. This 
was obtained using the basic setup but with a wave analyzer 
instead of the correlator. The frequency sweep was made very 
slow in order to get a good average. The main use of this 
figure was to provide the proper limits of integration for 
Equation (2.26). The limit of integration became a quite 
important parameter since it modifies greatly the form of the 
normalized beam pattern; thinking in terms of a transducer, 
it is known that higher frequencies make the beam narrower 
and the low frequencies have a wider beam, so looking at 
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Figure 3.6 Near field experimental setup. 
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Figure 3.2, the actual beam pattern, it is logical to assume 
that the main contributions are given by the low frequencies , 
and that is exactly what the Frequency Spectrum shows in 
Figure 3.7. The lower limit of integration is assumed to be 
1100 Hz and the upper limit is 5000 Hz, the contributions of 
the other frequencies was small; this also showed that the 
speaker has a low frequency response, although it covers the 
important band of the human hearing. 

Another critical factor in the Near-Field measurements 
was the correlation distance, since its effect in the width 
of the beam is great. The symmetry of the source caused some 
difficulties; it was found that all points located at the 
same distance from the axis of the horn had a very good cor- 
relation, regardless of the distance between them, but along 
the same angular position a distance of 12 cm was enough to 
diminish greatly the spatial autocorrelation function. This 
suggested that the correlation distance for this experiment 
is a vectorial parameter, but was used like a scalar along 
this work. It is useful to point out that in a real situa- 
tion with an asymmetric source driven by several sources with 
different frequency bands , like a submarine in which each 
section is affected for different set of frequencies and also 
the physical configuration changing in space, the autocorre- 
lation distance has to be easy to measure experimentally and 
its vectorial characteristic should disappear leading to a 
scalar parameter. Otherwise the pattern must be built up 
from several sources. 
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Along this experiment a correlation distance of 12 cm 
was used. 

C. RESULTS FROM NEAR-FIELD MEASUREMENT 

The first step in the data processing was to get the 
average mean square particle velocity. Using the data col- 
lected in the experiment and the formula given by Equation 
(3.1) it was possible to evaluate the power at each point of 
the active area, this power was averaged, the entire process 
was done using an IBM- 3 60 computer program in fortran IV 
language. The value obtained in this program represents the 
mean square of the particle velocity expressed by the 
Equation (2.22), where the value of the autocorrelation 
function is at zero spatial and time shift. 

The approximation used for that evaluation was, of course, 
an additional source of inaccuracy in the experiment, but 
should not affect too much the final result, because what is 
really affected is the intensity field. The normalized beam 
pattern is practically unaffected since the relative power 
between each angular position was the subject of interest. 

In the case of trying to obtain a radiation pattern, which 

. 2 

IS beyond the scope of this work, then Uo becomes a very 

important factor since it affects the Power Level of the 

signal. 

The solution of Equation (2.26) is the actual prediction 
of the Far Field Beam Pattern. Again the computer system 
was used for performing numerical integration, a library 
sub-routine called DQATR was used. It was necessary to use 
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double precision in order to obtain a reasonable accuracy. 

The sub-routine uses the trapezoidal rule in connection with 
Rombergs principles for evaluating the approximation of the 
integral value of the function. 

The resulting computed Far-Field is shown in Figure 3.8 
in comparison with the actual Far-Field. Note that the 
agreement between the predicted and the measured beam pattern 
is fair in shape but slightly out in values, the actual beam 
pattern shows a wider beam that the predicted one. Several 
factors can produce that difference, the correlation distance 
could be one of the main factors due to its specially vari- 
able characteristic in this experiment, as it was explained 
before. Also the limits of integration assumed an important 
role, strong high frequency components could produce that 
narrow beam. 

In order to study the variation of the function with 
respect to the limits of integration more calculations were 
made reducing the upper limit to 4250 Hz, all other parameters 
were kept constant. The result of this calculation is shown 
in Figure 3.9. This time a better agreement was obtained 
between the actual and the predicted beam patterns. A wider 
beam was obtained but gtill showed a small difference with 
the actual beam pattern. 

The effect of the correlation distance is shown in Figure 
3.10. The limits of integration were kept at 4250 Hz and 
1100 Hz, and the correlation distance was increased to 15 cm 
in the calculation. It can be seen that a small variation in 
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Predicted and measured beam pattern 
(UL = 5000 Hz, d = 12 cm, LL = 1100 Hz) 
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Figure 3.9 Predicted and measured beam pattern 

(UL = 4250 Hz, d = 12 cm, LL = 1100 Hz) 
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Figure 3.10 Predicted and measured beam pattern 

(UL = 4250 Hz, d = 15 cm, LL = 1100 Hz) 



this parameter affects greatly the predicted width of the 
beam. This Figure also shows that the beam pattern is more 
sensitive to the correlation distance than to the limits of 
■integration. 
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IV. CONCLUSION 



The possibility of predicting the far field beam pattern 
from randomly vibrating surfaces has been demonstrated. It 
is feasible to obtain the intensity field and the radiation 
pattern from near field measurements of distributed sources. 
Although the theoretical model developed in Section II was 
tested experimentally with satisfactory results, more work 
needs to be done in both areas . The amount by which the beam 
pattern depends on each parameter makes it difficult to state 
a solid argument for describing the range of validity of each 
one . 

The computed beam pattern was obtained solving equation 
(2.26) using a computer program for solving the integral by 
numerical methods. The limits of integration were chosen by 
means of the Frequency Spectrum, which was taken using a wave 
analyzer with automatic frequency sweep. The maximum range 
considered was from 1100 Hz to 5000 Hz. It was assumed that 
the frequencies producing the field were inside of this fre- 
quency band. 

The mean square of the particle velocity, which has 
energy information, was shown to be related to the autocorre- 
lation function. A correlator was used for that measurement 
and the power contained in the signal was easily calculated 
with equation (3.1). Note that in this case a correlator was 
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used but a true RMS voltmeter can also be used to measure 
the power of the signal. 

The correlation distance was also measured with the 
correlator, a pair of microphones were located in the active 
area and the two dimensional (spatial and time) correlation 
function was observed. A correlation distance of 12 cm was 
used for calculation, but the measured value was found to 
vary with location of the microphones . The calculated beam 
pattern is very sensitive to this parameter. Its effects 
are bigger than those due to the limits of integration. 

Different configuration of sources should be tested ex- 
perimentally. Sources should be excited for different fre- 
quency bands in such a way as to produce an active region 
with a well defined correlation distance. This will solve 
one of the main problems of this experiment. Also it is 
suggested that the frequency spectrum of the configuration 
should not necessarily be continuous, as in this form the 
arrangement will be closer to real cases. 

The theoretical model does not consider the spectrum 
level of the signal for the different frequencies. This 
should be tested more carefully in order to state a rule for 
the limits of integration. In the case of discontinuities 
in the Frequency Spectrum the model fails to give a straight- 
forward solution, although it looks reasonable to perform a 
summation of the contribution of each frequency band. It 
was not possible to verify this concept with the experimental 
setup . 
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Agreement was generally good between measured and com- 
puted patterns, especially when the upper limit of integra- 
tion was reduced to 4250 Hz. The actual beam pattern was 
slightly broader than the computed pattern in all cases. 

This is not surprising since the far field beam pattern de- 
pends on the frequency response of the source and the limits 
of integration were chosen in a rather empirical way, so the 
computed beam pattern could include a bigger frequency band 
than the real one and the width of the beam would be too 
narrow. 

The technique developed here might also be useful in 
measuring the radiated noise pattern of submarines. Some 
additional considerations should be taken into account, and 
the theoretical model should be validated by actual measure- 
ment in anechoic tanks . 
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APPENDIX A 



COMPUTER PROGRAM FOR EVALUATING THE FAR FIELD BEAM PATTERN 
USING THE THEORETICAL MODEL DEVELOPED IN SECTION II 



PROGRAM TO EVALUATE THE BEAM PATTERN OF A HORN 
SPEAKER DRIVED WITH RANDOM NOISE. 

FOR PROGRAMING PURPOSES THE FOLLOWING SYMBOLS HAVE BEEN 
USED: 

Cl = CHARACTERISTIC IMPEDANCE OF THE MEDIUM 
D = CORRELATION DISTANCE. 

BETA = ANGLE WITH RESPECT TO THE ACOUSTICAL AXIS. 

T = AVERAGE OF THE MEAN SQUARE OF THE PARTICLE 
VELOCITY IN THE ACTIVE REGION. 

XL = LOWER LIMIT OF INTEGRATION. 

WN=UPPER LIMIT OF INTEGRATION. 

A = AREA OF THE ACTIVE REGION. 

C = SPEED OF SOUND IN THE AIR. 



CALCULATE THE INTENSITY IN THE FAR FIELD 



IMPLICIT REAL*8 (F-H,Q-Z) 

REAL LABEL/ » » / 

REAL*8 title (12) 

REAL*8 PI,A,AUX(999) ,d,c,ci,cs 
DIMENSION P(73) ,BETA(73) 

EXTERNAL F 

COMMON D,C,RBETA,SN,PI 
PI = 3.1^115926 
D = 0.12 
EPS = 0.001 
WN = 2.*PI*5000. 

NDIM = 999 
A =0.5**2 
T = 0.062 
Cl = 415. 
c = 343. 

R = 5.8 

XL = 2.*PI *1100. 

51 = CI*A* T *(D**2) 

WP = 4250. *2. * PI 

52 = (R**2) *(C**2) *WP*((2.*PI)**1.5) 
DO 10 I = 1 ,37 

BETA = BETA(I)*PI/ I 80 . 

SN = (DSIN(RBETA) )**2 
CS = DCOS(RBETA) 

CALL DQATR(XL,WN,EPS,NDIM,F,Y,IER,AUX) 
P(I)=(S1*(0S**2)*Y) / (S2*((1.+0S)**2)) 
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BETA(74-I) = -BETA(I) 

P(74-I) = P(I) 

WRITE(6,100) Y,IER,BETA(I) ,P(I) 

10 CONTINUE 

100 PORMAT(/»Y =' ,E16.7,5X, ’IER=’ ,I5,5X, 'BETA=’ ,P7.2, ’INT= 
1 ',E16.7) 



CALCULATE THE NORMALIZED AMPLITUDE OP THE PRESSURE PAR 
PIELD BEAM PATTERN. 



Z = P(37) 

DO 20 1=1,73 

P(I)= ((P(I)/Z)»*0.5) 

WRITE(6,200) BETA(I),P(I) 

P(I) = P(D* 2. 

20 CONTINUE 

200 PORMATC/, ’BETA =» ,P7 . 2 , 13X, ' NORMALIZED AMPLITUDE =’,P 
17.5) 

read(5,50) title 

50 PORMAT(6a8) 

CALL DRAW( 73 , BETA, P, 0 , 0 , LABEL, TITLE ,20. ,0.5, 0,4, 0,2, 
19, 8,0, LAST) 

write(6,4oo) last 

400 PORMAT(/»LAST =',I2) 

STOP 

END 



PUNCTION P(W) 

IMPLICIT REAL*8 (A-P,0-Z) 

COMMON D,C,RBETA,SN,PI 
BIP =0.5 
BPP = 4250. 

ERP = 2. 

WP = BPP * PI * BRP 

P =(W**2)/((DEXP(BIP*((W*D/C)**2)*SN))*(DEXP(BIP* 
1((W/WP)**2)))) 

RETURN 

END 
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